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Error  Probability  Characteristics  for  Multiple  Alternative 
Communication  With  Diversity,  but  Without  Fading 

Iniroduciion 

The  M-ary  character  error  probability,  P^„  for  orthogonal  multiple  alicrnative 
communication  with  D-Cold  diversity,  was  evaluated  in  reference  1  for  a  wide  range 
of  parameter  values.  It  was  presumed  there  that  the  slow  signal  fading  on  the  D 
diversity  channels  was  independent  and  that  the  received  signal  energy  was 
distributed  exponentially  (Rayleigh  fading  of  received  signal  voltage). 

!'  is  of  interest  here  to  reconsider  the  system  performance  for  the  case  where, 
although  the  receiver  was  designed  for  D-fold  diversity,  there  is  in  fact  no  signal 
fading.  The  fractionalization  of  the  received  signal  into  the  diversity  chanttcis  theti 
results  in  poorer  performance  than  had  the  signal  been  confined  to  one  chatmcl  and 
prtK'csscd  coherently.  This  situation  can  arise  naturally  in  practice,  as  for  example, 
as  a  result  of  multipath  arrivals,  or  it  can  come  about  intentionally  in  the  system 
design.  It  can  also  be  of  interest  in  testing  a  multiple  alternative  hardware  design 
under  cottirollcd  laboratory  conditions,  where  the  time  or  cost  of  simulating  actual 
fading  conditions  is  excessive. 

The  basic  framework  and  background  for  the  communications  technique  con- 
idered  here  have  already  been  presented  in  reference  1  and  will  not  be  repeated,  for 
l  ie  sake  of  brevuy.  The  reader  is  presumed  to  be  familiar  with  the  time-bandwidth 
di, rat  ion  and  separation  constraints  listed  in  the  above  reference,  particularly  pages 
2-4  and  appendix  B. 
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Definitions  and  Technical  Results 


The  source  selects  one  of  M  equi-probable  symbols  and  transmits  a  sij;n;il  on  D 
diversity  channels  (time  and/or  frequency)  to  a  receiver.  1  he  D  channels  devoted  to 
each  of  the  M  alternatives  are  disjoint  (in  time  and/or  frequency)  with  each  other 
and  with  the  channels  for  the  other  alternatives  (see  reference  I).  The  receiver 
employs  matched  filtering  on  each  of  the  totality  of  MD  cells  ol'  interest.  F-dr  each 
signal  alternative,  the  corresponding  D  envelope-squared  matched  filter  outputs  arc 
sampled  appropriately  in  time  and  summed.  The  largest  of  these  M  decision 
variables  is  theti  declared  to  be  the  transmitted  signal  alternative. 

The  additive  noise  at  the  receiver  is  assumed  to  be  white  Ciaiissian  o\cr  the  total 
band  of  the  possible  received  signals,  with  a  (single-sided)  power  density  level  of  N,, 
watts/H/.  It  is  shown  in  appendix  A  that  the  M-ary  character  error  pnrbability  of 
this  cotninunications  technique  is  given  by 


P  =  1 
e 


(I) 


Here  “deflection”  statistic  d.  is  given  by 


d*"  E 

_T  _  _T  _  total  received  signal  energy 
2  "  "  noise  power  density  level 


(2) 


where  E-p  is  the  total  received  signal  energy  on  the  D  channels;  l„(x)  is  the  modified 
Bessel  function  of  order  n  and  argument  x;  auxiliary  function  B„(x)  is  defined  as 


'0 

and  (rcfeicnce  2,  equation  6.5.1 1 ) 


,(x)  ^  L  ^ 

‘  k=0 


(3) 


(4) 


is  the  exponential  power  series  through  the  term  x". 

It  is  interesting  and  worthwhile  to  observe  from  (I)  and  (2)  that  the  exact  frac- 
lionali/ation  of  the  total  received  signal  energy  E,  amongst  the  D  diversity  branches 
is  immaterial  in  so  far  as  system  performance  is  concerned.  The  available  signal 
energy  E ,  could  be  divided  equally  in  the  D  branches,  or  it  could  be  concentrated  in 
just  one  branch;  iti  either  ease,  error  probability  P^,  is  identictil  (assuming  M  and  1) 
are  unchanged). 

for  M  =  2,  expression  (1)  can  be  evaluated  in  closed  form.  Several  alternative 
representations  for  this  binary  error  probability  are  given  iti  appendix  B,  tilong  wiili 
an  efficient  program  for  this  special  case.  The  progrtim  for  cvtthiation  of  (1)  in 
general  is  given  in  appendix  C'. 
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(•raphical  Resulls 

For  a  particular  value  of  M,  the  character  error  probability  in  (1 )  is  a  lunciion 
of  the  number  of  diversity  channels,  D,  and  the  total  energy-density  ratio  (signal-to- 
noisc  ratio)  E,/N^,.  We  have  plotted  P^.  on  a  normal  probability  ordinate  versus 
Ef/N,,  on  a  dB  abscissa  (i.e.,  10  log  E^/N^,),  with  D  as  a  parameter,  iti  the  range  10  ' 
to  10  '  for  P^,.  Figures  1-1 1  correspond,  respectively,  to  M  =  2.  3,  4,  8,  16,  32,  64, 
128,  256,  512,  1024.  The  curve  for  D  =  1  in  each  figure  eorrespt)nds  to  the  NO 
FADING  results  iti  reference  I,  figures  1-7. 

We  hase  discovered  no  simple  rule-of-thumb  for  the  additional  signal  eneig\ 
required  to  maitttain  a  fixed  P^.  as  D  is  increased,  which  covers  the  whole  range  of  I) 
and  M.  Nor  did  vse  find  a  simple  expression  for  P^.  when  it  is  small  10  '),  despite 
considerable  effort.  The  best  we  can  do  is  to  observe  that 

Pg  «  (M  -  1)  P^2 

and  then  use  result  (B-IO)  for  the  binary  error  probability,  P^,,.  D  is  arbitrary  in  (5). 
Further  development  of  this  approach  appears  in  appendix  B. 
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Ej/No  in  dB 


Figure  11.  M-ary  Character  Error  Probability  for  M  =  1024 
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Correct  Decision  with  a  Threshold 


The  processing  system  above  presumed  that  one  of  the  M  signal  alternatives  was 
always  transmitted.  A  more  general  situation  occurs  when  the  additional  ease  of  no 
signal  present  at  all  is  allowed.  Then  the  largest  of  the  M  decision  variables  is 
compared  with  a  threshold.  There  are  three  types  of  probabilities  of  interest  for 
signal  present  (reference  3,  pages  2-5);  they  are  the  probability  of  a  tnissed  decision, 
the  probability  of  an  ineorreei  decision,  and  the  probability  of  a  eoneet  decision. 
The  last  is  given  by  a  slight  generali/at  ion  of  (I)  to  (reference  3,  equation  5) 


where  V  is  the  ihresluikl  value. 


f  (ir  signal  absent,  we  then  have  false  alarm  probability  (reference  3.  equation  2) 


FA  ■  ^  ’  Probjmax(yj.  ...  ,  <  V  for  d,j.  =  O] 


D-1 


=  1 


-  lo 


M 


=  1  -  [bp_j(vV2)]  .  (7) 

Mere  we  ttsed  the  limit  of  the  integrand  of  (6)  as  di-^O  atid  (3).  Result  (7)  agrees 
with  r.’ferenee  3,  equation  25,  as  it  must,  since  both  receivers  are  processing 
identical  noise-only  channel  outputs.  No  numerical  results  oti  this  more  geneitil 
situation,  (6)  attd  (7),  are  presetiied  here. 

tAtension  (o  Fading  Signal 


All  the  earlier  results  in  this  report  have  pertained  to  eotisiatit  amplitude  signal 
Ci)tnponents  on  the  D  diversity  branches.  Now  we  presume  tiuit  there  is  slow  signal 
fading,  meaning  that  the  total  received  signal  energy  F'  |  is  a  randotn  vartahle.  Itt 
particular,  we  let  ettergy-density  ratio  (signal-to-noise  ratio) 


K  -  — 
‘t  2 


have  probability  density  futtctioti 

,v 


R,j,  exp^-R.j,/p) 

r(v  -r  1) 


for 


(S) 


(9) 
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V^o 

V  +  1  ”  V  +  1 


,  V  >  -  1  . 


Here  Ej  is  the  average  received  total  signal  energy  on  the  p  branches,  and  v  is  a 
constant,  not  necessarily  an  integer. 

This  very  general  model  of  fading  was  presented  and  uiili/ed  in  reference  4,  pages 
1 1 -1 2,  in  the  investigation  of  a  maximum  likelihood  detector.  For  example,  the  four 
signal  fading  cases  considered  by  Swerling  (reference  5)  arc  subsumed  by  (9)  for 
particular  choices  of  v;  sec  table  1 . 

Table  I.  Particular  Cases  General  Fading 


V 

0  D-l  1 

213-1 

Swelling  Case 

1  2  3 

4 

Cases  I  and  3  correspond  to  common  signal  fading  on  a||  D  brandies,  whereas 
cases  2  and  4  correspond  to  independent  and  identically  distributed  signal  fading  on 
each  of  the  D  branches.  In  particular,  case  2,  v  =  D  -  1,  corresponds  to  ex- 
ponential  probability  density  function 

PjCW  ■•=  «P(--^)  tor  «  >  0  .  ^  .  1,11 

for  the  received  energy-density  ratio,  E,/N„,  on  one  branch:  this  is  Rayleigh  fading 
of  the  received  signal  voltage  amplitude  on  each  branch.  Equation  (9)  is  a  D-fold 
convolution  of  (11)  when  v  =  D  -  I.  Case  4,  v  =  2D  -  I,  corresponds  to 
probability  density  function 


p,(R)  =^R,xp(..^r) 


for  R  >  0  ,  Rj  “  » 


for  the  received  energy-density  ratio,  E,/N^,  on  one  branch.  Equation  (9)  is  a  D- 
fold  convolution  of  (12)  when  V  =  2D  -  1. 

The  average  character  error  probability  is  evaluated  in  appendix  D;  it  is  given  by 

p;  -  1  -  [(D-D!  (U  *  D"**] 

'  fo  rh-')  •  ,13) 

This  single  integral  can  be  evaluated  numerically  for  any  v  of  interest;  of  course,  M, 
D,  and  Ei/N„  need  to  be  specified  also  for  this  numerical  procedure. 

As  a  special  instance  of  fading  (9),  consider  case  2  in  more  detail: 
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V  =  D  -  1, 


Ilj.  E^/N^ 

M  =  -p-  =  — g —  =  average  energy-density  ratio  per  branch.  (14) 
Then,  if  we  use  reference  2,  equation  13.6.12,  and  (3)  above.  ( 1 3)  specializes  lo 


P  =  1  - 
e 


Jo  [Jo  '’“‘"’J 


(1^1 


where 


.D- 


= 


fo.  .  >  0  . 


(D  -  1):  (p  -r  1)‘ 


D-1  ,  , 

r„s  X  exp(-x)  ^  » 

=  •  (b  -  1) !  foi-  X  >  0  ,  (If,) 

are  ilie  proliabiliiy  densiiy  functions  of  the  decision  variables  formed  by  the  sum  i>l 
n  envelope-squared  filter  outputs  for  signal-present  and  noise-only,  respectively. 
Result  (15)  is  identical  to  reference  I,  equation  C-1.  It  is  what  we  wiuild  have  giMten 
had  we  averaged  the  signal  probability  density  function  in  (1)  —  the  function 
multiplying  the  bracket  —  with  respeet  to  the  signal  strength,  before  expressing  I’ , 
itt  integral  ftirnt.  Itt  fact,  a  more  general  integral  than  (15)  has  already  been  en- 
eotiniered  and  evahiaied  in  reference  3,  equation  26,  where  a  thresholding  operation 
was  included  itt  the  receiver  processor  for  the  fading  signal. 

For  M  -  2,  general  fading  lesult  (13)  can  be  evaluated  in  closed  form;  these 
special  cases  are  presented  itt  appendix  D,  especially  (D-5)  and  (D-8). 

Discussion  and  Summary 

I  his  report  has  addressed  the  problem  of  numerically  assessing  the  addiiiotial 
signal-to-tmise  ratio  (i.e.,  ettergy-density  ratio)  required  when  a  tionfaditig  sigtial  is 
brokett  ittto  D  (itneqital)  cotuponents  and  conibiticd  incoherently  .  Some  relaieil 
studies  into  the  cost  of  imperfections  of  receivers,  t>r  the  cost  of  lack  of  knowledge 
of  the  received  signal  structure,  are  presented  and  evaluated  in  references  6  and  7. 

Since  die  performance  of  the  processor  considered  here  for  the  nonfading  signal 
depends  only  on  the  total  received  signal  energy,  regardless  of  how  fractionalized  it 
may  be  among  the  diversity  channels,  these  results  also  apply  to  a  siiiiaiion  where 
imeeriainiy  in  the  lime  of  arrival  or  the  doppler  shift  of  the  received  signal  ettiiscs 
the  receiver  to  search  over  several  channels,  even  linmgh  the  particular  received 
signal  alternative  occupies  only  one  channel.  The  D-fokI  summation  lot  each  M-aiv 
alierniiiive  leads  Kt  a  kvss  in  perlormance;  stated  alternatively,  the  cost  nf  the  im- 
ceriainiy  is  additional  received  signal  energy  required  in  order  tiv  maintain  the  saute 
quality  of  performance.  1’^.. 

For  a  large  signal-to-noise  ratio,  the  M-ary  character  errvn  iiiob.tliiliiy  is  given 
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appioximaicly  by  (5)  in  terms  of  the  binary  error  probability,  P^.,.  The  usual  ap¬ 
proach  is  then  to  use  (A-3)  with  M  =  2,  or  to  use  (B-4),  to  evaltiale  However, 
the  situation  frequently  arises  w'here  the  characteristic  functions  of  the  dccisiini 
variables  can  be  evaluated  fairly  easily,  whereas  the  corresponding  probability 
density  functions  and/or  cumulative  distribution  functions  cannot.  In  that  case,  an 
alternative  representation  of  P^,,  directly  in  terms  of  the  pertinent  characteristic 
functions  would  be  useful.  This  problem  is  addressed  in  appendix  E.  with  the 
follow  ing  results  for  the  binary  error  probability; 

+ 

•  117, 

Here  f^U)  and  f|(4)  are  the  characteristic  functions  of  the  noise-only  and  signal- 
plus-noise  decision  variables,  respectively.  C_  is  a  contour  in  the  complex  4-plane 
along  the  real  axis,  with  a  small  downward  identation  at  4  =  0;  C+  is  a  similar 
contour  indented  upward  at  4  =  0.  The  last  form  in  (17)  can  be  particularly  useful 
for  numerical  evaluation  of  for  some  characteristic  functions. 

Another  alternative  to  evaluation  of  M-ary  error  probabilities  is  to  use  bounds 
which  are  simpler  ti)  compute  than  the  exact  result.  This  is  attractive  if  the  bt'unds 
;irc  tight,  at  least  for  the  range  of  signal-to-noise  ratios  and  error  probabilities  of 
interest.  Some  general  bounds  on  error  probability  (which  were  not  used  here)  arc 
presented  in  appendix  F  lot  a  situation  that  includes  interference  as  well  as  noise. 
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Appendix  A 

Derivation  of  M-ary  Character  Error  Probability 


To  avoid  duplication  of  effort,  wc  will  rely  heavily  on  reference  S.  pages  2,  3,  8,  9 
and  appendixes  A  and  B,  especially  (B-7)-(B-12).  We  obtain  directly  the  probability 
density  function  of  the  signal  decision  variable  (sum  of  D  filter  outputs)  as* 


and  that  for  any  one  of  the  M-1  noise-only  decision  variables  as 

=  i  ^''p  (- 1)  y  >  0  •  <A-2) 


A  correct  decision  is  yielded  only  if  x  is  greater  than  all  M-l  indepettdeni  ni'ise 
variables;  the  probability  of  this  event  is 


P 


c 


dx  pj(x) 


dy 


PoW 


M-l 


(A-3) 


But  the  integral  on  y  in  (A-3)  and  (A-2)  is  given  by  (3)  as  B,)  |(x/2).  Then  bv  use  of 
(  A-l )  and  the  substitution  X  --  y-,  (A-3)  becomes 


P 

c 


exp 


(- 


4/2) 


.D-1 


/■ 

•'0 


J 

dy  y 


-y^/2 


M-l 


(A-4) 


1  he  M-ary  character  error  prr'bability  is  P^.  =  I  -  P^.. 

Several  checks  on  (A-4)  arc  available.  First,  for  the  trivial  case  M  =  1.  (A-4) 
yields  P  =  1,  which  agrees  with  the  fact  that  there  are  then  no  noise-only  channels 
to  worry  about  at  all.  Second,  as  d  ,“*■0,  (A-4)  reduces  to  P^  -  I  /M;  this  agrees  u  ith 
the  observation  that  we  have  a  eomplclcly  random  selecliim  with  zero  signal-ti>- 
noise  ratio.  I'liird,  for  D  1 .  (A-4)  agrees  with  reference  1 .  equation  9A  (when  the 
latter  is  corrected  for  the  typographical  omission  of  the  factor  exp(-[-.  j/N^,))- 

A  program  lor  the  evaluai ion  of  P  ,  =  1  -  P^,  via  (A-4)  is  given  in  appendix  C. 


*A  ‘‘defleeiion”  interpretation  of  is  available  from  reference  8,  equations  (B-IO),  (B-7). 
(B-4),  (A-4),  and(B-5). 


A  1  A  2 
Reset  SC  Blattk 
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Appendix  B 

Binary  Krror  Probability 


This  appendix  will  deal  exelusively  with  ihe  special  ease  (d  M  2.  a  binaiA 
decision.  Then  (I )  becomes,  with  the  ideniiricaiit>ns  in  (3)  and  (4). 


Tt)  evaliiaie  the  integral,  we  used  reference  9.  equation  6.631  I:  the  final  irtins- 
formaiion  to  (H-l)  employed  reference  2,  equation  13.1.27.  An  alicrnaiise  ex¬ 
pression  to  (B-1 )  is  available  via  reference  2,  equation  13.6.9: 


P 


e2 


1 


D-1 
kaO  2^ 


(D-1) 

k 


(B-2) 


If  we  expand  the  polynomial  in  d^/4  that  oceurs  in  (B-1 )  and  (B-2),  tind  interchange 
sums,  we  obtain 

exp(-d^/4)  D-1  fdJ/4)"  D-l-n  (D  +  n)  . 

p  =  — L  .1  f  Y  y.  —■ - ^  . 

®  2®  2"  n!  2^  j:  (B  3) 

All  three  expressioiis  above  lor  are  finite  sums  of  positive  quantities  and  arc 
reasonable  frtr  computer  evahtaiir)n,  even  for  fairly  large  1). 

Another  expression  for  P^-,  =  I  -  P  ,  is  available  by  recalling  (A-l)  attd  (A-2) 
and  changing  (A-3)  to 

-‘>'PoW  I  ■‘■‘P.W  ■ 

But  from  (A-l),  upon  making  the  substitution  x  1-  attd  using  icfcicttcc  10, 
cqimtion  I,  we  fittd  for  the  inner  iniegral  in(B-4), 
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D-1 


*2  .2' 
t  . 


_/•*"  I  t  + 


)-i(V) 


(B-5) 


Then  (B-4)  yields 


c2 


/: 


D-l 


a,  t.:  (  ^1/2) 

2  (D  -  1):  ^  ' 


L 


=  [2°-^  (D  -  1).']  du  exp(-uV2)  Qp(d^.  u)  . 

•\i  this  pciini,  we  will  use  (he  new  integral  resiili 
dx  exp(^-p^x^/2)  Qw(a,  bx) 


2'^’^(N  -  1): 
2N 


,M+N-1 


1  -  - 


-A 


N-l/M+N-r 


2 

D  +  b 


®N-l-n 


w  here 


A  =  ^ 


2  2.2 

p  +  b 


(B-7) 


(B-8) 


1  his  inlegral  and  several  other  new  integrals  of  Q^,  fimelions,  along  with  their 
derivations,  will  be  presented  fairly  soon  by  this  author  in  another  NUSC  teehnieal 
I  eport . 

I  he  use  of  (B-7)  and  (B-8)  in  (B-6)  yields 
exp 


e2 


1  ("■')  ('M 


(B-9) 


I  vpanding  the  partial  exponential  expansion  in  powers  of  d^/4  aeeording  to  (4), 
and  interehanging  summations,  we  obtain  for  the  binary  error  probability 


e2 


exp(-d^/4)  D^l  (d2/4)  D-^-k  pD-1 

20-1  ~  “  iT*  ”  “  \  n 

2  k=0  n=0  ' 


(B-IO) 


I  bis  result  is  very  similar  to  (B-3),  and  of  course  yields  ideniie;il  numerical  values; 
however,  a  direct  verification  of  the  transformation  that  would  take  the  inner  sum 
iH  (B-IO)  to  the  inner  sum  in  (B-3),  or  vice  versa,  has  not  been  discovered. 
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(B-10)  is  a  very  efficient  form  for  numerical  evaluation,  because  the  summand  ot 
the  inner  sum  (on  n)  does  not  depend  tm  the  index  (k)  of  the  outer  sum;  this  ad¬ 
vantageous  property  is  not  true  of  (B-3).  Thus  for  a  given  D,  we  can  compute  the 
partial  sums  on  n  in  (B-10)  by  starting  at  k  =  I)  -  1  and  merely  adding  one  ad¬ 
ditional  term  to  the  inner  sum  for  each  downward  unit  step  in  k;  these  D  partitti 
sums  are  ail  stored  and  then  used  for  the  outer  sum  starling  at  k  =  0.  I  lms,  whttt 
appears  to  be  a  double  sum  in  (B-10)  actually  can  be  computed  by  means  of  two 
single  sums  of  D  positive  terms.  Furthermore,  the  innei  sum  on  n  in  (B-IO)  can  be 
done  once  for  a  value  of  D  of  interest,  and  then  used  repeatedly  lor  dificretn  values 
of  d^.  A  program  incorporating  these  observations  is  listed  ;it  the  end  of  this  ap¬ 
pendix.  All  of  the  results  thus  far  apply  to  arbitrary  values  of  energs -dciisii y  ratio 
Et^N,,. 

One  reason  for  dwelling  on  the  evaluation  of  the  binary  error  prrrb;ibilii>  is  tluit 
lor  large  signal-lo-noise  ratio,  i.e.,  large  energy-density  ratio  E  |/N  the  M-ars 
character  error  probability  is  given  approximately  by 

P  J5S  (M  -  1)  P  -  for  small  P  .  (B-11) 

e  e2  e 

I  bis  result  may  be  seen  to  be  valid  when  it  is  reali/cd  that  the  probabilitv  of  ;m\ 
particular  noise-alone  decision  variable  exceeding  the  signal  decision  rariahlc  is  ,i 
rare  cvcttl  when  P^,  is  small.  Thus  the  probability  of  two  in  more  noise  variables 
exceeding  the  signal  variable  is  very  rare  indeed,  and  we  need  only  account  foi  the 
single  noise  variable  case.  Since  these  events  are  disjoint  and  there  are  M-1  such 
variables,  (B-ll )  follows  immediately. 

For  large  d^/2  =  an  asymptotic  form  for  P^.,  can  be  developed  from 

(B-10);  express 

exp(-d^/4)  (d^/A) 

^e2  ‘  (D  -  1) : 


2D(D  -  1)  ^  (D  -  1)  (D  -  2)  (2D"  -  D  -r  1)  ^ 
2  /  -,  \  2 


dj/A 


Therefore,  asymptotically. 


(di/^y 


(B-12) 


exp^-d^/4j  {d^/4') 


D-1 


'’e2  (D-1):  as  d,j,  -  “  . 


(B-1.3) 


But  reference  to  (B-12)  reveals  that  in  order  for  (B-13)  to  be  a  good  approximation 
to  P..,,  rather  than  just  an  asymptotic  result,  the  requirement  d,  '>>  1)  must  be 
satisfied.  This  requirement,  however,  is  typically  too  large  a  signal-io-noise  latio 
requirement  to  be  practically  useful.  We  have  lound  it  necessary  to  resort  to  (B-10) 
for  the  binary  error  probability  for  the  range  of  values  of  itiieresi  bet  e. 
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A  program  (or  the  binary  error  probability  P^.,  in  BASIC  Cor  the  Hewlcii- 
Packard  9845  is  presented  below.  It  utilizes  (B-10)  et  seq.,  as  diseiissed  above.  The 
variables  are  self-explanatory  and  consistent  with  the  earlier  text. 


;  '  E 1  r,  "  Pe  tor  D-fotd  divers;  t;.',  but  oo  fadirij.  TP  e.473 

tC  Ii=12  !  ri>0  NUMBER  OP  IlIvERSITY  CHHNNEL'i 

iO  EtriO  =  20  I  Et,-,0>0  TOTRL  SIGNftL  ENERGY  SIHGLE-illiEIi  NOIEE  DENSITi' 

30  r'l=D-l 

40  D2  =  I(»i; 

50  f3  =  Erri0  2 

Eo  DIM  i';i00.i 

70  REDIM  SsDi:' 

S0  ■3'::Di:i=T=l 

■30  FOR  .T  =  l  TO  Dl 
100  T  =  T*  C  Ii2-  J  >  J 

110  El  di-  j>=S':d-J)+t 

120  f(E;-:T  J 

1  3  0  S  =  Si '  0  ;.'i 

140  T=1 

150  FOR  K=1  TO  Dl 
1  r'  0  T  -  T  ♦  R .  K 

170  S  =  StTRS<K'' 

ISO  ne;;t  k 

r30  Ps  =EXP  ( -R  ■' *'3  •  2  1,  D2-l> 

200  PRINT  'll  = :  0 , Et  IT'S  = " ; Et ri0,  "Re  «'’;Pe 
210  END 


=  20 


Pe  =  3.  4513133'3i37E-03 


E »  fi  0 
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Appendix  C 

Program  for  M-ary  Character  Krror  Probability 


The  expression  for  P^.  is  given  by  (1).  Since  we  cannoi  integraie  lo  0°,  we  lerminaie 
ihe  integral  at  y  =  L.  The  error  incurred  in  doing  this  is  given  by 


exp 


f  dy  y°  e-y  lD-l(‘^Ty) 


< 


dy  y 


^D-l(‘^Ty) 


exp[-  1(l  -  d^)^] 

\V  (l  -  d^)  ■ 

We  used:  the  lact  that  B|,(x)  is  upper-bounded  by  I,  as  may  be  seeti  froiti  (3);  the 
asymptotic  behavior  of  l|,(x)  b'l"  large  x  as  given  by  reference  2,  equation  9.7. 1 ;  and 
an  integration  by  parts  lo  yield  the  final  form.  For  given  values  of  D  and  d,  (M  has 
dropped  out  of  this  error  bound),  we  slen  !  up  by  units  of  1  from  the  value  d,  4  I 
until  (( -1 )  is  sufficiently  small. 

We  chose  error  10  in  the  enclosed  program.  The  reason  that  the  integral  in  (I ) 
must  be  evaluated  very  accurately  is  that  we  must  subtract  the  integral  result  Irom  1 
in  order  to  gel  the  error  probability.  If  we  want  to  evaluate  P^,  accurately  in  the 
range  of  10  ^  then  we  need  (he  integral  to  8  or  9  decimal  places  to  ciumteraci  the 
effect  of  subtraction  from  1  required  by  form(l). 


Cl 
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Wc  numerically  integrate  (1)  for  y  in  the  interval  (0,  1.)  via  the  Trapezoidal  rule, 
with  automatic  halving  of  the  interval,  until  a  stable  result  occurs.  It  is  shown  in 
reference  1 1,  pages  55  and  58,  that  the  Trapezoidal  rule  is  excellent  for  integrands 
w  hose  odd  derivatives  at  the  end  points  of  integration  are  equal.  Kor  the  application 
here,  the  integrand  of  (1)  may  be  shown  to  behave  as  y-i>M  as  y-*0+  .  Thus  a  high 
number  of  odd  derivatives  are  zero  at  y  =  0+  and  virtually  zero  at  y  =  L  for  large 
DM,  meaning  that  the  correction  terms  to  the  Trapezoidal  rule  vanish.  In  fact,  (1) 
was  also  integrated  via  Simpson’s  rule  and  via  Simpson’s  rule  with  end  correcliott 
(reference  12,  pages  414-418);  neither  performed  as  well  as  the  Trapezoidal  rule, 
w  hich  is  incorporated  in  the  program  for  P^.  below.  Because  of  the  very  small  values 
lor  the  exponentials  and  the  very  large  values  for  the  Bessel  function,  it  was  found 
necessary  to  first  evaluate  the  logarithm  of  the  integrand  of  (1)  and  then  ex¬ 
ponentiate  it  prior  to  its  use  in  the  Trapezoidal  rule. 
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J&O  IiEF  FNS'  V  '  f  IHTEGfttiND 

31  u  COM  D,  D1  ,  Ml  ,  Dt  ,  Et  ri0 

320  T  =  .  '.  ♦  V  *  V 

330  E=1-EXP‘'-T  '•^FMECT.m  > 

340  IF  B'=0  THEN  RETURN  0 

350  ■5  =  Ii+-L0G  ‘  V  Dr  )  -  T  +  FNI  nx  1  og<  Dt  *Y,  Dl  >+M  1  ♦lOG<  B  >  -  Et  no 

3t0  RETURN  EXP^S'^Dr 

370  FNEND 

330  I 

330  DEF  FNEwl.N)  !  SUM  OF  X'K'K!  FROM  K  =  O  TO  N 

400  ■3  =  T  =  1 

410  FOR  K=1  TO  N 

420  T=T^X  1 

430  5=3+T 

440  NEXT  P 

450  RETURN  S 

400  FNEND 

470  ! 

430  DEF  F  N I  n  ;■  1  o  g  <  X ,  N  )  •  L  0  G  U  F  I  n  ^  x  > 

430  IF  ;,>0  THEN  520 

500  IF  N=0  THEN  RETURN  1 

10  IF  N;0  THEN  RETURN  0 

20  Fl=.5*X 

30  F  =  1 

40  FOR  1=2  TO  N 

50  F=F*I 

tO  NEXT  I 

70  f=n»log‘:fi''-logcf;‘ 

30  F(  =  Fl*fl 

30  S=T=lE-87 

oOO  FOR  1=1  TO  500 

S10  T=T*F|/CI*CN+I)^ 

S20  S=S+T 

630  IF  T*2Ell;=S  THEN  EGO 
t'4  0  next  I 

€•50  OUTPUT  0;"50O  TERMS  AT  ";X;N 
€•€0  I  n>.  1  .:.g=F+LOG S +200.  32430  309 
€70  RETURN  Inxlog 

€80  FNEND 
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Appendix  1) 

Average  Error  Probability  for  Fading  Signal 

The  energy-density  ratio  variable  was  defined  in  (8),  and  its  probability 
density  was  given  in  (9)  for  a  fading  signal.  From  (I),  the  error  probability  expressed 
in  terms  of  R  is  given  by 

exp(-llj.) 

D-1 

K) ' 

(D-l) 


■fhe  average  error  probability  is 


■  /„"  ““t  pCV)  ".W  • 

Substitution  of  (9)  and  (D-l )  in  (D-2)  and  interehange  of  integrals  lead  to 


(D-2) 


P  =  1 
e 


•'0  .  u  rcv  +  1) 

K) 


(D-3) 


In  the  inner  integral,  let  Rj  =  I-/2  and  use  reference  9,  equation  6.631  I.  Sim¬ 
plification  leads  to 

1-1 

-1, 


P  =  1 
e 


r  ixV+ll  r"  .  2D-1  -7^2 

-  1^2  (D  -  1) :  (y  +  1)  J  j  dy  y  e  ^  f 


(D-4) 


Use  of  the  variable  t  =  y’/2  in  (D-4)  yields  the  result  already  quoted  in  (13).  Results 
(D-4)  and  (13)  hold  for  general  v  and  M.  Wheti  v  is  chosen  as  the  special  value  D-l, 
the  resulting  average  error  probability  is  given  by  ( 15)  and  (16). 

Specialization  to  IVI  =  2 


Now  we  let  parameter  v  in  (9)  be  general,  but  we  set  M  =  2.  Thus  we  w  ill  develop 
the  average  binary  error  probability  for  a  gen.;ral  signal  fading  model.  Upon  itse  o' 


D-l 
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(3)  and  (4),  (13)  yields  for  M  =  2, 

-,-1 


■  [»  -  1) 


!  (U  *  !)'’*'] 


(P.  *  V  (ii-L'f  f(d  .  1  .  D  .  k;'  D;  ^ 

Cw  .  2)“  fro  ^  Ap  •  2/'  ^  ’  -  •  P  * 


D-l-\)  D-1 


(U  +  2) 


1  liis  finiic  sum  of  Ciaussian  hypergeomciric  fiincdons  is  not  loo  useful  in  general. 
However,  for  ilie  special  ease  v  =  D  -  1  considered  in  (14)-(16),  (D-5)  immediately 
yields 

—  1  /D-l+k\/M  +  ^ 

P  T  =  - fT  V  (  |( - =  )  for  V  =  D  -  1  ; 

(U  *  2)°  kto  W  A'"  " 


this  last  result  can  also  be  derived  directly  from  (15)  and  (16)  lot  M  -  2.  by 
teletetice  to  (3)  and  (4). 

A  much  more  appealing  result  than  (D-5)  for  average  binary  error  probability  P  , 
for  general  v  is  attained  if  we  start  with  the  binary  no-fading  result  iti  (B-9)  and  use 
idetttificatioti  (8): 


„  exp  (.R,,/2)  D^l  ^20.1^  D-^-n  ^ 

~  "  22D-1  n^O  V  "  /  h  ■  <’ 

Then  the  average  binary  error  probability  is,  by  use  of  (D-7)  and  (9), 

^e2  ^  Jo  p('*t)  **62 

L  J  n=0^  "  '  k=0  k!  2*^ 

•  /o  “P  (-  T  ■  ■ 

r,2D-l  /.  V.\'"‘T'  l'2D.l\  “-fr"  ‘P  *  ')k  /  p/2  ) 

■p  r?)  J  „?oUJ  k?o  (r-rsTj) 


U  ■  2 


D-2 
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This  closed  form  result  is  a  finite  sum  of  positive  quantities  and  holds  for  general 
values  of  u(>-l).  The  summations  are  very  similar  to  those  encountered  in  non¬ 
fading  binary  error  probability  (B-10).  A  program  for  the  evaluation  of  (l)-8), 
which  lakes  advantage  of  the  observations  made  under  (B-IO),  is  given  below .  I  or 
V  =  1)  -  1,  (D-8)  must  reduce  to  (D-6);  wc  have  not  discovered  the  iratisformaiion 
that  accomplishes  this  corroboration,  but  we  have  confirmed  the  cqu;ilii> 
numerically  for  v  =  D  -  1,  as  well  as  for  the  alicrnaiive,  more  general  rcsuli  (l)-5). 
No  numerical  results  for  (D-8)  have  been  presented  here. 


1  '  fi'csfigc  biriiiM  Pe  for  D-fold  diMsrsity  with  fading. 

10  D=l£  !  D.:ei  HUMBER  OF  D IVERS I T'-  CHMNHE.L:. 

^■0  !  Ho>-l  MEftSURE  OF  FFtHIHG,  EQS.  9-10 

50  I  Mo>  =  0  MEftSURE  OF  S  IGHftL-TO-HO  I SE  RftTIO 

40  PRINT  "D  =";r,“Hu  =";Nu,“Mu  =“;Mo 

50  Dl=r-l 

60  Ii2- D*2 

70  ft-Mu  s2  +  I'K4) 

SO  DIM  S<100> 

90  PEDIM  S<Dl> 

100  SsDri=T=l 

110  FOR  J=1  TO  D1 
120  T  =  T*'.D2- J:-/ J 
150  SsDl- J)=S<:D- J>  +  T 

140  NEXT  J 

150  S  =  S<.0> 

160  T=1 

170  FOR  K*1  TO  D1 
ISO  T  =  T  *  (N  u  +  K  '  *  ft  K 
190  s=3+t*S(t;;j 

200  NEXT  K 

210  Pe-^S  '  2'  '  D2  -  1  >  ■* '•  1  + .  5*Mu > * < +  1  >  > 

220  PRINT  "Pe  =";Pe 
250  END 


D  =12  Nu  =  6.1  Mu  =  5.7 

F  s  =  6  .  1  0  6  1  5  5  I  6  S  6  5  E  -  0  3 


l)-3/D-4 
Reverse  Blank 
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Appendix  F 

Binary  Error  Probability  Via  Characteristic  Functions 


Let  p„(y)  and  P|(x)  be  the  general  probability  density  liineiions  ol  a  noise-only 
decision  variable  and  a  signal-plus-noise  decision  variable,  respectively.  1  he  binary 
probability  of  error  is  then 


“  Prob(x  <  y) 


dx  pj(x) 


dy  p^(y) 


(L-1) 


Now  by  relerenee  13,  equation  6,  the  exceedance  probability  in  (L-l)  is  given  in 
terms  ol'theeharaeieristic  fiiiieiion  f, (4) according  to 


dy  PgCy)  =  1  -  Pq(x) 


1  J_ 

2  *  i2ir 


^  ^  f  (O  expC-i^x) 


=  ife  jT  -f  ^0^^^  expC-iW  .  (E-2) 

The  integral  with  a  slash  on  it  is  a  principal  value  integral;  the  contour  C_  is  a 
contour  on  the  real  axis  of  the  complex  {-plane  except  for  a  small  downward 
identation  at  {  =  0.  Employment  of  (E-2)  in  (E-1)  yields 

dx  Pj(x)  ik  T 


■  TH  /  f  'o'O  ^I'-O 


’■m  I  fV-«  h'«  .  (i-j, 

where  C  ,  is  indented  upward  at  4=0.  Both  contours  C' ,  start  at  4  =  -<»  and  go 
to4  =  + 

Expression  (E-3)  gives  binary  error  probability  directly  in  terms  of  charac¬ 
teristic  functions  f^fl)  and  f|(4),  by  means  of  a  contour  integral.  Since  characteristic 
functions  satisfy  the  relation 

f(-C)  =  f*a)  for  real  £.  .  (1-4) 


E-1 
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an  altcrnaiivc  real  integral  form  for  (E-3)  is 


^  -f  00 

I  f 


E-2 


(-Q} 
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Appendix  F 

Bounds  on  Error  Probabilily 

We  consider  the  situation  where  there  are  three  different  variables  upon  which  a 
decision  as  to  the  transmitted  multiple  signal  alternative  must  be  tnade.  1  lies  are, 
respectively,  the  signal,  noise,  and  interference  variables.  In  pariieulai ,  lei  random 
decision  variable  S  denote  the  output  of  the  signal  channel;  S  ma\  iisell  be  the  sum 
of  several  diversity  branch  outputs.  Let  N  be  the  output  of  the  nuiximifiii  o(  the 
noise  channels  outputs,  of  which  there  may  be  many;  each  ni>ise  ehannel  onipni 
may  itself  be  the  sum  of  several  diversity  branch  ouipnis.  And  let  I  be  the  onipm  ol 
the  maximuni  of  the  interference-plus-noise  ehannel  rniipuis,  i>f  which  there  may  be 
many;  each  inierferenee  ehannel  output  may  itself  he  the  sum  of  seseial  diversiix 
branch  outputs. 

I'or  independent  signal,  noise,  and  inierferenee  ehannel  imipnis,  the  probabiliis 
of  a  eon  eel  ileeision  is 

=  Probd  <  S,  N  <  S)  =  y*dx  p^(x)  P.(x)  P^^Cx)  ,  (I  D 

where  p,(x)  is  the  probability  density  function  of  random  variable  a,  ami  L  ,(v)  is  its 
eorresponding  cumulative  distribution  function.  W'e  also  express 

P.(x)  =  Probd  <  X)  =  1  -  Q.(x)  , 

Pj^(x)  =  Prob(N  <  x)  =  1  -  (}^(x)  ,  (1-2) 

where  Q,,(x)  is  the  exceedance  probability  of  random  variable  a. 

The  probability  of  character  error  is 

’’e  =  ^  ■  ‘’c  = 

=  y*dx  Pg(x)[Q^(x)  +  {J^(x)  -  Q^(x)  Q„(x)] 

=  Probd  >  S)  +  Prob(N  >  S)  -  Probfl  >  S,  N  >  S) 

(1-3) 

I  he  last  quantity  in  (F-3)  is  generally  quite  small  for  practical  ctiscs  of  imcicsi.  Hut 
in  any  event,  we  have  the  upper  bound 

Pg  <  Probd  >  S)  +  Prob(N  >  S)  (and  Pg  <  1)  •  (1  -f) 


In  addition,  since 
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